J ftp suppression in the threshold model and QGP formation time 
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In the QGP motivated threshold model, in addition to the normal nuclear absorption, J/tp's 
are subjected to an additional "anomalous" suppression. We have analysed the recently published 
PHENIX data on the participant number dependence of the nuclear modification factor for J/ip's 
in Au+Au collisions and extracted the anomalous suppression required to explain the data. At mid 
rapidity J ftp's are anomalously suppressed only above a threshold density n c =3.73 fm -2 . The for- 
ward rapidity data on the otherhand require that J ftp's are continuously "anomalously" suppressed. 
The analysis strongly indicate that in mid rapidity J ftp's are suppressed in a deconfined medium. 
Using the PHENIX data on the participant number dependence of the Bjorken energy density, we 
have also estimated the QGP formation time. For critical temperature T c =192 MeV, estimated 
QGP formation time ranges between 0.06-0.08 fm/c. 
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Recently, PHENIX collaboration have published their 
measurements of the centrality dependence of J ftp sup- 
pression in Au+Au collisions at RHIC energy, ^/i=200 
GeV Data are taken at mid rapidity (\y\ < .35) and 
at forward rapidity (1.2 < \y\ < 2.2). In most central 
Au+Au collisions, J ftp's are more suppressed at forward 
rapidity than at mid rapidity. Suppression factor is ~ 3 
at mid rapidity and ~ 6 at forward rapidity. 

There is growing consensus that in central Au+Au col- 
lisions at RHIC, a deconfined state of quarks and gluons 
(QGP) is produced. It is expected that a deconfined 
medium, if produced in Au+Au collisions will leave its 
imprint in J ftp production. Long back, Matsui and Satz 
3 predicted that in a deconfined medium, binding of a 
cc pair into a J ftp meson will be hindered, leading to 
the so called J ftp suppression in heavy ion collisions Q . 
However, J ftp's are also suppressed in a nuclear medium. 
Inelastic interactions of J ftp's with the nucleons can dis- 
sociate them. Suppressed J ftp production not necessarily 
imply a deconfined matter formation. At RHIC energy, 
it has been further argued that rather than suppression, 
J ftp production will be enhanced 0, H[. Due to large 
initial energy, large number of cc pairs will be produced 
in initial hard scatterings. Recombination of cc can oc- 
cur, enhancing the charmonium production. However, as 
mentioned earlier, PHENIX data do not show any indi- 
cation of J ftp enhancement. 

PHENIX data on the centrality dependence of J ftp 
suppression has been analysed in several models, e.g. e.g. 
comover model _@, statistical coalescence model @, the 
kinetic model 0, S| or the QCD based nuclear absorp- 
tion model [l|. None of the models give satisfactory de- 
scription of the experimental data. Recently, we have 
analysed [l(| the PHENIX data in the threshold model 
and found that it do explain the PHENIX mid rapidity 
data, but not the forward rapidity data. In the present 



* E-mail: akc@veccal.ernet. in 



paper, we refine the analysis and extracted the threshold 
density required to fit the mid rapidity and the forward 
rapidity PHENIX data. Extracted threshold density is 
then used to obtain a physical parameter, the QGP for- 
mation time. We find that QGP formation time is quite 
small, r 0.06 - 0.08/m. 

Blaizot et al [ill, |TJ] , proposed the threshold model to 
explain the NA50 data on anomalous J ftp s upp ression in 
158 AGeV Pb+Pb collisions at SPS energy [jj] . Thresh- 
old model tries to mimic the sudden melting of J ftp in 
a deconfined medium. In the model fate of a J ftp de- 
pend on the local energy density. If the energy density 
exceeds a critical value, the inter quark potential can not 
bind a cc pair into a J ftp. It is also assumed that "lo- 
cal" energy density is proportional to 'local' transverse 
density. Then if the " local" transverse density exceeds a 
critical or threshold value, deconfined matter is formed 
and all the J ftp's are completely destroyed (anomalous 
suppression). One must remember that the anomalous 
suppression is in addition to the "conventional nuclear 
absorption" . The model neglects the transverse expan- 
sion of the system. It is implicitly assumed that J ftp's 
are absorbed before the transverse expansion sets in. 

In the threshold model, number of J ftp mesons, pro- 
duced in a AA collision, at impact parameter b can be 
written as, 

N%*( h) = N W J ^ sTl // (s)T e// (b _ g) 

x S anom (b, s) , (1) 
where T e ff(b) is the effective nuclear thickness, 



T eff {b) = dzp(b,z)exp{~<7 abs dz/p(b,z/)), 

J —oo J z 

(2) 

dabs being the J/-0-nucleon absorption cross-section. For 
the density we use the Woods-Saxon form. 
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= TT FT pv v / rf3r ^ r ) = A ( 3 ) 

1 + exp((r — R)/a) J 

with i?=6.38 fm and a=0.535 fm. 

S anom (b,s) in EqQ]is the anomalous suppression fac- 
tor introduced by Blaizot et al. G3 • Assuming that 
all the J/ifj's get suppressed above a threshold density 
(n c ), the anomalous suppression can be written as, 

S ano m{b, s) = 6(n c - n(b, s)) (4) 

where n c is the critical or the threshold density. n(b, s) is 
the local transverse density. It was observed [llj that by 
smearing the threshold density by a small amount, one 
can obtain better fit to the data, but at the expense of 
an additional parameter (A), 



However, experiments do indicate otherwise. In Au+Au 
collisions, J/ip's are more suppressed at forward rapidity 
than at mid rapidity. In the threshold model, such a de- 
pendence can only be accommodated if parameters of the 
model, (Tabs, n c and A depend on the rapidity variable. 
We thus separately fit the mid rapidity and the forward 
rapidity PHENIX data to extract those parameters. Be- 
fore we proceed further, we would like to note that it is 
not unnatural to have rapidity dependence on the critical 
parameter n c . For example, it is well known that the crit- 
ical temperature of the confinement-deconfinement phase 
transition depend on the baryon density of the system. 
Mid-rapidity region is essentially baryon free while at for- 
ward rapidity baryon content is non-negiligible. Rapidity 
dependence of the critical parameter n c will then implic- 
itly account for the baryon dependence of the critical 
parameter. 



S a r. 



,(b,s) = - [1 -tanhA(n(b,s) - n c )] (5) 



Critical ingredient of the threshold model is the " local" 
transverse density. At impact parameter b and at the 
transverse position s, local transverse density it can be 
obtained as, 



n(b,s) = T A (s)[l-exp(-a NN T B (s-b))} 

+T B (b-s)[l-exp(-a NN T A {s))] (6) 
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FIG. 1: PHENIX data on the centrality dependence of J/ip 
suppression in Au+Au collisions, at mid rapidity and forward 
rapidity. Best fit to the data in the threshold model is also 
shown. 

With anomalous suppression defined as in Eq[5l the 
threshold model have three parameters, J/V'-nucleon ab- 
sorption cross-section a a0S1 the threshold density n c and 
its smearing A. In the threshold model, J/ip suppres- 
sion do not depend explicitly on the rapidity variable. 



FIG. 2: Anomalous suppression as extracted from the 
PHENIX data on the centrality dependence of J/ip suppres- 
sion. The solid and the dashed lines corresponds to mid ra- 
pidity and the forward rapidity data. 

In Fig[Tl PHENIX data are shown. Data points are 
few and it is not judicious to fit all the three parameters 
simultaneously. With the Glauber model of nuclear ab- 
sorption, we first fit the few peripheral collision (up to 
N par t=150) in mid rapidity data and extract the J/ip- 
nucleon absorption cross-section, a a b s — 4.39 ± 0.7imb. 
The value is larger than the estimated u a bs— 1-3 mb 
fl4| in d+Au collisions. In peripheral Au+Au collisions, 
J/^'s are more suppressed than in d+Au collisions. With 
a aos fixed from peripheral collisions, we fit the full data 
set to find the threshold density n c and its smearing 
A. Best fit to the mid rapidity data is obtained with 
n c = 3.73 ± 0.29/m- 2 and A = 8.96 ± 9.72/m 2 . The 
solid line in Fig Q] shows the fit. The quality of fit is very 
good. For the forward rapidity data sets also we use the 
Co6s=4.39 mb. As seen in Fig[TJ at extreme peripheral 
collisions, J/ip suppression in mid and forward rapidity 
is similar. Best fit to the forward rapidity data set is ob- 
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tained with n c = 2.96 ± 0.42/m" 2 and A = .99±.94/m 2 . 
Here again, as shown in Fig[T] the quality of fit is very 
good. While the threshold model do explain the cen- 
trality dependence of J/ip suppression at mid rapidity as 
well as at forward rapidity, the anomalous suppression 
{Sanom) required for the two data sets are widely differ- 
ent. In Figl^l we have shown the anomalous suppression 
Sanom as required by the mid and the forward rapidity 
data. At mid rapidity, true to the spirit of the threshold 
model, anomalous suppression shows a step like behavior. 
At mid rapidity, J/ip are " anomalously" suppressed only 
above the threshold transverse density rt c =3.73 fm~ 2 . 
But at forward rapidity J/t/j's are continuously "anoma- 
lously" suppressed. Even though the model fits the data, 
the spirit of the model is lost. Step function like anoma- 
lous suppression in mid rapidity give strong indication 
that at mid rapidity, J ftp's are suppressed in a decon- 
fined medium. Continuous anomalous suppression at for- 
ward rapidity on the other hand indicate that J/ip's are 
possibly suppressed due a mechanism not related to the 
confinement-deconfined phase transition. However at for- 
ward rapidity, J '/ip suppression is more complex than en- 
visaged in a simple Glauber like model. Simple Glauber 
model can not explain cither of the two data sets. 
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FIG. 3: Solid circles are the PHENIX data on the partic- 
ipant number dependence of Bjorken energy density times 
QGP formation time (tebj) [13] ■ The curve is the maximum 
transverse density that can be achieved in Au+Au collisions. 



Before we proceed further, we would like to note that 
the threshold density as determined here represents the 
upper limit. Threshold model neglects some very im- 
portant effects, e.g. (i) feed back from ip' and \ states 
and (ii) transverse expansion. A considerable fraction 
(~ 40%) of J/ip's are from decay of %p' and \ states (l5j . 
That part is completely neglected here. Threshold den- 
sity for anomalous suppression of higher states, ip' and 
X should be less than that for a J/ip. Then presently 
estimated threshold density n c represent an upper limit. 
Additionally, at RHIC, model studies indicate that in 



the deconfined phase, the system undergoes significant 
transverse expansion [iH ]. The local transverse density 
is a key ingredient to the Threshold model. In an ex- 
panding system, local transverse density will be diluted. 
J ftp's, which are anomalously suppressed in a static sys- 
tem, may survive in an expanding system due to dilu- 
tion. Then, the presently estimated threshold density 
will again represent an upper limit. 

We now try to connect the estimated threshold den- 
sity with some physical parameters like threshold en- 
ergy density or temperature above which J/ ip are anoma- 
lously suppressed. As mentioned earlier, threshold den- 
sity is assumed to be proportional to energy density. 
If the proportionality factor is known, we can estimate 
the threshold energy density above which the J ftp's are 
anomalously suppressed. As given in Eq|5l local trans- 
verse density is a function of the impact parameter (b) 
and the transverse position (s). For collisions between 
two identical nucleus at impact parameter b, maximum 
transverse density is achieved at the transverse position 
s = b/2. In Fig|3]we have plotted the transverse den- 
sity n™ ax (h) — n p (h,s = b/2) as a function of partici- 
pant number. n™ aa; increases with the collision centrality. 
If in a collision with participant number N part , decon- 
fined matter is produced and J ftp's are anomalously sup- 
pressed, at the minimum n p nax should exceed the thresh- 
old density. As seen from Figj3l estimated threshold den- 
sity, n c =3.73 ± 0.29 fm~ 2 corresponds to Au+Au colli- 
sions with participant number N part — 199.6igg'2. 

Experimentally one estimate the initially produced en- 
ergy density by measuring the total transverse energy 
Et and using an estimate for the initial reaction volume. 
In the Bjorken model with longitudinal boost-invariance, 
the energy density is obtained as, 



£bj 



1 dE T 
tAt dy 



(7) 



where t is the formation time, At is the overlap area and 
dEx/dy transverse energy per unit rapidity. QGP for- 
mation time is an important parameter. Experimental 
determination of energy density then depends strongly 
on the estimate of the initial time. PHENIX collabo- 
ration have measured the transverse energy Et- Since 
QGP formation time is not known, they have tabulated 
the Bjorken energy density times the formation time as 
a function of the participant number. In Fig[3l we have 
shown the PHENIX data on the participant number de- 
pendence of the tebj [17j. Like n™ aa: , te bj increases 
as the collision centrality increases. PHENIX data in- 
dicate that a collision with participant number N part = 



3.98It'^GeV/m- 



199.6^552, corresponds to te th 
Eth is the threshold energy density above which J ftp's 
are anomalously suppressed. Corresponding threshold 
temperature (Tth) can be easily obtained using the re- 
lation, e = gQCP^jT 4 with (7qgp=47.5, for a QGP with 
three flavors. 
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TABLE I: Threshold temperature (Tth) above which 
J/ip's get dissociated and critical temperature (T c ) for the 
confinement-deconfinement transition for various choices of 
QGP formation time (r). 



T 


Tth 


T c 


(fm) 


(MeV) 


(MeV) 


0.02 


558.6±^j 


266.01^ 


0.04 


469.7Ig 1 '5 


223.71^ 


0.06 


424 4+ 24 ' 9 

^^•^-46.6 


202.1^22.2 


0.08 


395.0-ttj 


188.1^20.6 


0.10 


373.6+Hl 


177.911°* 


0.12 


356.91™ 


170.011°;? 


0.14 


343.4I3?; 2 


163.5l? 7 6 9 


0.16 


332 1+ 19 - 5 

"•J^- 1 — 36.5 


158.2l? 7 3 4 


0.18 


000 c + 
O^^.U 35 4 


153.6l? 6 ° 9 


0.20 


314.11^ 


149.6I? 6 8 4 



Tth is the temperature above which J/tp get dissoci- 
ated. Lattice based potential models indicate that in a 
deconfined medium, at the critical temperature T C1 in- 
terquark potential is not sufficiently screened to dissoci- 
ate J/ip's- Model calculations indicate that J ftp's can 
survive upto a temperature of 2.1T C In tablelH for a 
choice of formation time r, we have tabulated the thresh- 
old temperature (Tth) and the critical temperature (T c ). 
For formation time varying between 0.02 fm to 0.2 fm, the 
critical temperature varies from 150 MeV to 265 MeV. 
Critical temperature for the confinement-deconfinement 
transition has been accurately estimated in recent lattice 
calculations, T c ~ 192±7±4 MeV [II]. As seen from ta- 



ble [fl it corresponds to formation time r ranging between 
0.06-0.08 fm. The time is considerably smaller than the 
thermal equilibration time r eq 0.6 fm [16]. QGP is 
produced early in the collisions. 

To summarise, in the QGP motivated threshold model, 
we have analyzed the PHENIX data on the centrality de- 
pendence of J ftp suppression in Au+Au collisions. In the 
threshold model, in addition to the normal nuclear ab- 
sorption, J/'i/j's are anomalously suppressed, such that, if 
the local transverse density exceeds a threshold density 
n c , all the J ftp's are absorbed. In a careful analysis, we 
have extracted the threshold density required to explain 
the mid rapidity and the forward rapidity PHENIX data. 
Mid rapidity data are well explained in the model with 
threshold density n c = 3.73 ± 0.29/to~ 2 . The data re- 
quire very small smearing of the threshold density, A = 
8.96 ± 9.72/m. The forward rapidity data on the other 
hand require very large smearing, n c = 2.963 ±0.42/m -2 
and A = 0.99±0.94/m. Very large smearing required for 
the forward rapidity data defeat the essence of the thresh- 
old model which tries to mimic the sudden onset of J / tp in 
a deconfined medium. We conclude that J/ip suppression 
at forward rapidity, though more complex than envisaged 
in the Glauber model of nuclear absorption, do not indi- 
cate a deconfinement phase transition. J '/ip suppression 
at mid rapidity which require sudden on set of anoma- 
lous suppression above the threshold value n c =3.73 ± 
0.29 /to -2 , possibly indicate a deconfined matter forma- 
tion. Using the PHENIX data on participant number 
dependence of Bjorken energy density times the forma- 
tion time, we have estimated the QGP formation time as 
t pa 0.06 — 0.08/m for critical temperature T c 192 
MeV. 
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